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ABSTRACT 


The mechanics of the oblique impingement of an axisymmetric 
jet have been investigated. The kinematic constraint of no penetration 
requires the jet to spread and/or to deflect vertically in order to provide 
a sufficient cross sectional area for the mass flux. The selection between 
these two effects is established by the dynamical effects represented in 
the three-dimensional momentum equation. The velocity field, as 
expressed by the convective acceleration of the mean flow is primarily 
balanced by the three-dimensional pressure distribution above the 
plate and tailored by the Reynolds stress distribution. It is shown that 
an alternative description of the pertinent phenomena may be constructed 
by considering the vorticity transport equation. The gradients of the 
(easily measured) two-dimensional surface pressure field are source 
terms for the lateral and longitudinal vorticity components and the 
azimuthal vorticity structure in the approaching jet is reoriented to 
create a deformation in the isotach pattern which is characteristic of 
secondary flows. 

The jet responds by deflecting vertically for a sufficiently slow 
approach to the plate; for this, the intersection of the conical free jet 
isotach with an imaginary reference plane characterizes the interaction. 
The spreading effect is important for larger angles and/or closer 
spacings. The stagnation point is located near the upstream edge of the 
interaction and is physically displaced from the maximum pressure 
point on the plate. 
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1 . INTRODUCTION 

The mechanics of the oblique impingement of an axisymmetric jet 
on a plane surface are examined in detail in this report. Figure 1 shows 
a schematic drawing of the problem under consideration and the coordinate 
system used to describe the flow field. The kinematic features of the 
flow above the plate are examined in the context of the conservation of 
mass, the vorticity of the jet, and the vorticity introduced by the jet- 
plate interaction. The dynamic features of the flow are examined in terms 
of the surface pressure distribution and the cause- effect relationships 
which exist between the pressure and velocity /vorticity distributions. 

This report represents the primary results of the second year's activities 
supported by the NASA grant NGR 23-004-068. 

This study was motivated by the externally blown flap (e.b.f. ), 
a configuration under consideration as a means of gaining large lift 
coefficients for STOL aircraft. For this application problem the airfoil 
extends no more than five or six jet diameters beyond the jet; hence, the 
extensive data for the inclination angle a and the height above the plate 
h/d are restricted to the limited streamwise domain of 0 ^ x/d ^ 6 for 
the static pressures. With the exception of the pressure data for a = 30, 
and a = 60 degrees, the present investigation is relevant to the flow 
resulting from the interaction of the propulsion jet with the main airfoil. 

The information herein is appropriate to an ( over- or unde* - wing configura- 
tion. A schematic of the latter is shown in Figure 1 The documented 
flow is an approximation of that which would be presented to the wing- 
flap juncture. The 30 and 60 degree conditions are representative of the 
direct interaction of the propulsion jet with the flap. 

The data base used to establish the mechanics of jet-plate inter- 
action is comprised of (1) surface static pressures measured at 0. 1 inch 
increments in y at discrete x locations, (2) velocity measurements from 
0. 1 inch increments in y at discrete z positions and at integer x/d 
locations from 0 to 5. These data have been used to create quantitative 
measures of the important phenomena occurring in the flow. Analytical 
considerations have been used to develop these measures and to infer 
significant relationships or implications which demonstrate the nature of 
the phenomena occurring in this flow field. 

The data base which serves as the primary resource for this 
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report was used in the First Annual Report to define measures of the 
flow field which were relevant to the externally blown flap problem. 

The measures and analyses of the present report are to define the basic 
character of the flow. This is considered to be complementary to the 
information of [l]. The identification of the oblique jet mechanics will 
contribute to several aspects of the externally blown flap development 
activity. Since the oblique jet impingement is the simplest form of the 
e.b.f. flow field, the basic aspects of the two flows will be similar. The 
mechanics identified in the present study cam serve as the basis for the 
interpretation of the e.b.f. flow field response to the complicating factors 
of (1) an external streaming flow, (2) the curvature and a finit : extent of 
the impingement plate (the airfoil), and (3) the presence of the deflecting 
flap. The mechanics identified herein will also allow an inference of the 
minimum spacing for multiple jets such that their flow fields show no sig- 
nificant interaction. The turbulent motions w'n*ch are responsible for the 
acoustic emissions from the impingement region are related to the mean 
flow field characteristics identified herein. Subsequent studies of the 
turbulence structure will be guided by these results. 

A general review of the literature pertinent to this problem was 
presented in [l] and no additional studies directly relevant to this problem 
have been identified. Consequently, a literature review will not be included 
herein. A brief description of the experimental facility is given in 
Section 2. Section 3 presents a comprehensive statement of the mechanics 
of the mean flow in the obliquely impinging axisymmetric jet. Specifically, 
Section 3. 2 identifies the spreading and curvature effects; Section 3. 3 
considers the vorticity aspects of the motion; Section 3.4 identifies a 
generalized framework for the development of a computing scheme; 

Section 3. 5 introduces the additional phenomena of the stagnation 
point and the y- component vorticity considerations; and Section 3.6. 
reiterates the important characteristics identified in the First Annual 
Report. Section 4 presents the detailed information which can be extracted 
from a consideration of the centerline pressure values. Section 5 presents 
a detailed consideration of the stagnation point: specifically. Section 5. 2 
presents the analytical considerations which demonstrate the consequences 
of separate physical locations for the stagnation and the maximum pressure 
points; Section 5. 3, the experimental results; and Section 5.4, the 
relationship of the isotachs and the stagnation point. The appendices are 
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used to communicate the pertinent results available in [l]. 


2. EXPERIMENTAL FACILITY 

The large volume of data considered in this report (and in [l ]) 
is essentially made possible by the on-line IBM 1800 digital computer 
facility. The quantitative evaluation of the several measures of the 
experimental data, i. e. integrations and the preparation of level curves 
(isotachs and isobars) is made feasible by the availability of the hot-wire 
and pressure transducer data on punched card This data processing 
facility and the experimental flow system are described in detail in [1]. 

Two nozzle configurations were used to examine the influence of 
the initial conditions of the jet-plate interaction. A fully developed flow 
was developed in a 2 inch I. D. , 12 ft long tube; this is referred to as the 
fully developed (f. d. ) condition. A second nozzle was used to approximate 
a uniform (u) flow. A nominally flat-top velocity profile was created 
with this nozzle and these results will be referred to as the uniform (u) 
condition. Details of this nozzle are given in Figure 2 along with an 
exit velocity profile. 

The velocity data which were used to create the isotach patterns 
(selected isotachs have been reproduced in Appendix A for reference 
purposes) were obtained with a single vertical wire (parallel with z). 
Consequently, these readings are most appropriately referred to as 
approximately (i. e. discounting pitch effects) the magnitude of the 
velocity in the x-y plane of the wire. Since the x- component is dominant 
except at the edge of the jet where u(x,y, z) << u(x,o, z), the reasonable 
approximation is that the hot-wire reading is the x-component of velocity. 
This distinction is necessary when the various flux integrals are formed 
or when the comparison with analytical flux predictions, e. g. the Reichardt 
model are examined. See Sections 4. 3 and 4.4 of [l], respectively, for 
these processed forms of the velocity data. 

A special series of experiments was performed to investigate 
the stagnation point. These experiments involved the flow, traverse, 
and data processing systems described in [1]. The appropriate details 
of the test procedure will be covered along with the presentation of the 
experiments. 
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3. THE MAJOR PHENOMENA OF THE OBLIQUE IMPINGEMENT OF AN 
AXISYMMETRIC JET 

3.1. Problem Statement 

The purpose of this section is to consider the dominant phenomena 
which occur in the jet-plate interaction. Subsequent sections will consider 
appropriate details of the separate phenomena. 

A useful conceptualization of the jet- plate interaction is based 
upon the no- penetration constraint imposed the plate on the motion of 
the jet fluid. That is, there is no z component of velocity in the plane of 
the plate. In its absence, the jet fluid which would have penetrated through 
the plane of the plate will pass a given x location by spreading laterally 
along the surh.ce of the plate and by increasing the area available for mass 
flux by causing the jet to curve upward. In the context of the kinematics 
of the problem, it is only necessary that the appropriate mass flux occurs 
across any x = constant location; the relative importance of spreading 
and curvature is dependent upon the dynamics of the jet-plate inter- 
action. That is, the pressure distribution within the jet determines the 
curvature and spreading effects. It is not possible to identify a cause 
and effect relationship based upon the conservation of mass equation 
since this is not the governing phenomena for the problem. 

The above considerations allow a somewhat more precise 
statement of the desired description of jet-plate interaction mechanics. 
Specifically, it should identify the magnitude of (1) the jet spreading and 
(2) the curvature of upward deflection effects; and identify the factors 
which are responsible for the selection process which results in the 
balance between these two effects. Clearly, the surface pressure distri- 
bution is mechanistically related to these two aspects of the motion in the 
jet. However, the surface pressure distribution itself is also an impor- 
tant aspect of the interaction process. Somewhat less obvious is the 
relationship between the surface pressure distribution and the x- 
component vorticity flux through an x = constant plane. This voriticity 
component is an important factor in the boundary layer control 
characteristics provided by the flow which bleeds between the lower surface 
of the airfoil and the flaps. The relationship between the pressure 
and the x- component vorticity will be examined in detail. The 
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curvature and spreading measures will first be examined; a model of 
the jet-plate interaction will then be presented. 

3. 2. Spreading and Curvature Effects in an Obliquely Impinging Jet 
3.2.1. Jet spreading 

From conservation of mass consideration, it is argued in Section 
3. 1 that the jet will both spread laterally along the surface of the plate, 
and curve upward. It is the purpose of this section to quantitatively 
evaluate these two effects. 

In order to evaluate the spreading effects, some knowledge or 
reliable estimate of a reference or intersection contour of the jet on the 
plate is necessary. Such a reference contour can be established by the 
intersection of the free (uniform or fully developed) jet with an imaginary 
plane in the same relative position as the plate. More specifically, if 
one considers the axisymmetric surface defined by the condition u(x, r) = 
constant, the intersection of this surface with the plane of the plate defines 
a contour. From the data of the present study, it is possible to quite 
accurately determine the functional relationship for the average r(x) 
values such that u(r,x)/u(o) = 0. 1, 0.3, 0.5, 0. 7, 0. 9; the average is 
based upon the large ensemble of experimental conditions with the following 
characteristics: (1) different a and h/d values, (2) an undisturbed 
isotach pattern, and (3) r values which are themselves smoothed over 
multiple y traverses. The 0. 1 isotach (u/u(o) = 0. 1) was chosen for this 
purpose since it best describes the outer boundary of the jet. Only the 
cases which demonstrated conservation of momentum were used to 
establish the radii of the 0. 1 isotach for the five downstream x/d locations. 
Table 1 presents these r(x) values. The linearly increasing values define 
a cone. The intersection of a cone with a plane defines an ellipse, 
parabola, or hyperbola depending upon the relative magnitudes of the 
cone and the intersection angles. These generated contours are termed 
free isotach intersections. The isobar and the free isotach intersection 
contours are shown in Figures 3 to 27. The shape of the two curves is 
remarkably similar which allows several measures of the interaction to 
be established. These are discussed below. 

The zero isobar is considered to be a significant measure of the 
outer jet boundary. That is, the zero isobar can be interpreted as the 
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limiting extent of the dynamic effects in the flow. The small negative 
pressures beyond this isobar are a result of the entrainment effects. It 
is reasonable to assume, but it has not been shown, that the zero isobar 
is related to the stagnation line where the jet fluid, flowing laterally 
outward along the plate, meets the entrained fluid. The stagnation pressure 
of the entrained fluid would be approximately (clearly no greater than) the 
atmospheric value. Since the free isotach intersection characterizes 
the size of the undisturbed jet in the plane of the plate, the ratio of this 
contour width to that of the zero isobar can be used to describe the 
relative jet width. If the isobar family of curves were parametric, then 
the x location used to measure the relative width would not influence the 
magnitude of the ratio since the free isotach intersection is defined by 
a cone and a plane. The isobars, especially the zero isobar, are 
apparently not members of a single family; therefore, there is some 
influence of x- location on the relati\ dth measure. The arbitrary 

decision was made to evaluate the width measurements at the x location 
of the maximum pressure. Since the zero free isotach grows at a 
different rate than the 0. 1 free isotach, the distance from the appropriate 
point at the nozzle exit to the maximum pressure location is considered 
to be pertinent supplementary information. Table 2 presents the percent 
width and nozzle-to-plate data taken from the isobar contours. 

The relatively constant percentage width values in Table 2 for the 
three to nine or twelve degree cases are considered to be quite significant. 
Contrary to what might be reasonably expected, it appears that there is 
very little systematic spreading of the jet with increasing q values. 
Considering a nominal value of 0. 7 for the ratio of these t.vo jet width 
measures, the majority of the data for both nozzles is within ± 0. 05 or 
= i 7 percent. Within this domain of cases, there are some second order 
systematic and non- systematic variations. The systematic deviations are 
dependent upon the increasing a values and probably the distance to tne 
maximum pressure location. The scattered or non- systematic variations 
may be caused by geometric, flow, or measurement condition abnormali- 
ties (for example, an improper a value in the experiment or a zero drift 
in the pressure transducers). The significant deviations from this 
cluster about 0. 7 are considered to be the u: a = 12, h/d = 1. 5; a = 15, 
h/d = 1; a = 30, h/d = 2,3; and a = 60, h/d = 2 cases in which a 
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pronounced and systematic trend is observed where the width of the jet 
is considerably larger than the free isotach intersection. These 
differences in the width ratios are apparently related to the much different 
interaction phenomena. Specifically, based upon these ratios and based 
upon the available isotach information, the cases with a nominal width 
ratio of 0. 7 appear to preserve the basic axisymmetric character of the 
jet and to confine the interaction phenomena to the region below the 
center of the jet. For the larger angle cases, the axisymmetric character 
of the jet is destroyed and the lateral flow along the surface of the plate 
is greatly enhanced as is the region of dynamic importance represented 
by the isobars. 

Table 2 also presents the magnitude of the normalized pressure 
of the isobar which aligns with the 0. 1 free isotach intersection contour. 

The pattern in which these values are essentially constant from three 
to twelve degrees is again apparent. A nominal value of 0. 27 and second 
order variations are indicated by the table. The percentage variations be- 
tween different cases in the magnitude of the isobar which aligns with the 
free isotach are larger than the percentage variations for the width 
measures. Also, the variations in the uniform cases are greater than in 
the fully developed cases and it can be observed that a pattern of larger- 
pressure for smaller-distance-to-impact is established, These are 
considered to be mechanistically related to the relatively steep velocity 
gradient at the edge of the uniform exit condition jet. The tabulated 
values are considered to be accurate within approximately ± 10 percent; 
the correspondence of free isotach/isobar contour was not exact and 
the equivalent isobar value often had to be obtained by interpolation. 

Large deviations from the nominal 0. 27 value are observed for 
the uniform cases of a =15, 30, and 60. The significantly larger pressure 
magnitudes for these cases and the much smaller percentage width 
measures noted above both signify qualitatively different interaction 
phenomena compared to the smaller angles. These differences, then, 
are to be related to the strength of the interaction. The strong interaction 
destroys the axisymmetric character of the upper portion of the jet; a 
weak interaction does not. 

A quite significant quantitative aspect of the isobar/free isotach 
alignment is that the factor of sin a , used to normalize the pressures, 
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allows the identification of a single nominal value over the range of 
a = 3 to 15 degrees. This is approximately a five-fold variation in the 
magnitude of sin a and hence a five- fold variation in the magnitude of 
the surface pressure. This simple dependence of the pressure on the 
jet angle would suggest that the z- component of the momentum flux is 
responsible for the elevated pressure in this region; this is in contrast 
to the centerline pressures which will show an additional effect associated 
with the overall jet and the local curvature effects. 

The isobars represent a measure of the dynamic width of the jet. 

The isobar field could not be wider than the jet flow but the converse 
could occur. That is, the inertia of the lateral motion induced by the 
pressure gradient within the active region of the jet-plate interaction could 
move the jet fluid beyond the zero isobar. In this event, the isobar would 
under- estimate the jet width. Figures 28 to 30 have been prepared to 
compare the isotach contours with the lateral pressure distribution for 
the uniform cases: a = 3, h/d = 0.75; a = 9, h/d = 1; and a = 15, h/d = 1. 
These figures demonstrate that the jet does spread beyond the zero isobar; 
however, the velocity beyond this boundary is relatively small. The 
above considerations are therefore restricted to the notion of the dynamic 
width of the jet. 

3.2.2. Curvature of the jet 

For the cases of weak jet-plate interaction, there are only minor 
changes in the spread as a function of a and h/d. Consequently, it is 
concluded that any significant differences between cases will be manifest 
in the curvature of the jet. The curvature of the entire jet can be computed 
from the surface pressure values. The assumptions and analytical 
considerations which make this possible are identified in the following; 
the curvature magnitudes will be examined following the analysis. 

Special experimental tests and an estimate of the maximum 
plausible influence of the static pressure was used in [1] to demonstrate 
that the x-component momentum flux was essentially constant for all 
cases. As a consequence of this, control volume momentum and 
moment- of- momentum analyses may be used to define the height above the 
plate z m (x) of the center of momentum flux J x (— J(o) cos a = constant). 

The quantity z m (x) can be considered to be the elevation of the jet (as if 
it were concentrated along a line). The curvature of the jet is given in 
terms of the curvature of z m (x) by the relationship 
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The expression for z m (x) developed in [1] will be used to evaluate the 
curvature. The expression from [l ] was made incomplete by the exclusion 
of a term which accounted for an entrainment effect. The corrected 
analysis is presented in Appendix B and the corrected form is 
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where is the entrainment velocity at the control surface and the 
control surface is defined in Figure 31 . The bracketed term is shown 
to be small for x/d i 5 in Appendix B. The derivatives for the 
curvature are evaluated as 
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Consequently 
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Although the curvature involves the second derivative of z m> it can be 
evaluated in terms of quantities which are quite accurately known, as 
shown by (5). 

An approximate relationship involving K is pertinent to the 
interpretation of the curvature data; specifically, since the maximum 
magnitude of dz m /dx is -tan a , the largest coefficient mul'iplying the 
numerator of K is [l+tan 2 a ] . ^ 2 Since 0. 9 £ (1 +tan 2 a )" 1 for OS 

Z 3 /2 

a £15 degrees, the variation of the denominator is neglected and [l+tan aj* / *1 
is considered to represent a satisfactory assumption. Consequently, 


K 3£ d 2 z /dx 2 
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( 6 ) 


The area under the K(x) curve is then 



(7) 


Since dz m (®)/dx = 0, the total area under the K(x) curve is approximately 
tan a . This may be used as a reference condition to compare the K(x) 
curves for a given a as a function of h/d. 

A second general observation may be made concerning the K(x) 
distributions. Unless the plate influences the jet at the e>.it of the nozzle, 
the curvature will vary from zero to a negative value (if the entrainment 
eLect is important) to positive values (for a sufficient duration to allow 
the area under the K curve to be = tan a ) and return to zero for large 
downstream distances. Th< centerline pressure will have the same 
general character; however, the pressure may be influenced by the local 
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effects in the flow whereas the curvature is a measure of the overall 
behavior of the jet. 

The curvature and the centerline pressure are presented together 
in Figures 32 to 45. The comparison of these two distributions is con- 
sidered to offer an effective measure of the curvature effects. If the 

two curves have the same character, it is concluded that the local 
curvature effects are not very important and that the jet is essentially 

turned en masse. The extent to which the two curves have a different 
character is used as an indication of the presence of local curvature or 
a relatively rapid jet approach to the plate. A shift in the x/d location has 
been applied to these (and subsequent plots) to align the 0. 1 free isotach 
intersection with the x/d = 0 location. This allows a better visual compar- 
ison between cases. The Ax/d values to accomplish the shift are 
presented in Table 2. 

The most important characteristic of these data is their slope as 
a function of the streamwise distance. The slopes were not evaluated 
since the usual errors associated with differentiation of data would have 
been encountered; the reader is asked to mentally compare the slopes 
of the K and p (x, 0,0) curves. The trend indicated by the uniform, 
a = 3 degree cases (h/d = 0.75 and 1) is characteristic of all the data. 

That is, the p(x,0,0) distribution reaches a maximum at a smaller x/d 
than does K(x) and this difference betv-een the two curves is less pro- 
nounced as h/d increases. This trend is quite pronounced for the fully 
developed, a = 3 degree cases. The condition wherein the centerline 
pressure reaches a maximum value before the curvature, reflects the 
dependence of p(x, 0,0) on the local flow condition. As h/d is decreased, 
the velocity gradients of the flow approaching the plate are increased and 
the surface pressure required to turn the flow locally is increased. The 
quantitative character of this observation is demonstrated by the plots. 

Figures 44 to 55 present a comparison of the curvatures for 
various h/d values at a given a . These plots are prepared such that 
the 0. 1 free isotach would be shifted to x = 0. The purpose of this 
comparison is to demonstrate the accommodation of the curvature to the 
different geometries. That is, the jet curves more rapidly for smaller 
h/d values. 

3.2.3. Summary of spreading/curvature effect 

The original question addressed in this section can now be answered. 
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It was noted that the constraints imposed by the conservation of mass 
required that the jet spread laterally along the plate and/or be curved 
upward because of the no penetration condition imposed by the plate. 

The relative spreading measures indicate that for the relatively small 
angles, a ? 12 degrees, the lateral spread of the jet is a second order 
effect. For these same cases, the curvature is a strongly varying function 
of a and h/d; that is, the jet primarily responds to the presence of the 
plate by turning upward. 

The larger angle cases appear to be qualitatively different. The 
upper portion of the jet does not retain its axisymmetric character and 
the relative width of the zero isobar to the 0. 1 free isotach is increased 
with respect to the smaller angle cases. The pressure distributions show 
a more peaked character with the relatively large maxima apparently 
accounting for the strong turning of the fluid at the upstream edge of the 
jet. 

3.3. Vorticity Considerations 

The upward deflection of the jet can be expressed in terms of the 
deviation of the z- component velocities in the jet from the equivalent 
values which would be realized if the plate were not present. Similarly, 
the spreading of the jet is associated with the excess y-component veloci- 
ties. The detailed velocity distributions are expected to be governed by 
the three-dimensional pressure distribution above the plate and tailored 
by the Reynolds shear stress distribution. The consideration of the jet- 
plate interaction mechanics based upon tl ese aspects iB reasonable in 
the context of the nature of the flow but such a description is quite 
difficult to relate to the extant or anticipate, experimental data; i. e. , 
there is a lack of p(x, y, z) data. An enuivalent, but conceptually quite 
different approach is to consider the flow in terms of the vorticity prccesses 
occurring in the jet-plate interaction. This has the implicit benefit that 
the important x-component vorticity will be directly involved in the descrip- 
tion of the problem. It has the major benefit that the three-dimensional 
pressure distribution in the jet flow does not enter into the description 
of the problem; rather, the surface pressure distribution, which is accurately 
known, enters the description as a source of vorticity. These considerations 
require the control volume description of the vorticity equations; these 
relationships are developed below. 

The Navier-Stokes equations for an incompressible flow can be 
operated upon to create the vorticity transport equation; specifically, 
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Performing a volume integral of this differential equation for a motion 
which ie steady in its time mean quantities and applying the Gauss 
Theorem to transform volume integrals of divergence terms into surface 
integrals yields (note that (“) denotes a time averaged value) 
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The term on the left is the net flux of the vorticity from the control volume 
by the convective action of the velocity field, jj,®n = ^ j£,»n + ^ ji°n 
where and & are the turbulent fluctuations about the mean). The first 
term on the right hand side of equation 9 is the production of vorticity by 
the amplification effect of stretching or reorienting vortex filaments. 

The third term can be expected to be important only at a solid surface 
where the viscous terms are important in an otherwise turbulent trans- 
port dominated flow field. The viscous term may be expressed in terms 
of the surface pressure gradients; this is developed in the following. 


d u 3 u) u 9 u 

FIT = 1 Tz + + k ITT 


( 10 ) 


Since u = (dx/dy - 3 v/fl z) and u = (9 u/9 z - 9v/9x), the vorticity 
x y 

terms at the surface of the plate may be written as 


<*> x (x, y, 0) 
Wy(x, y, 0) 



z-0 



z-0 


(11a) 

(lib) 


These quantities may, in turn, be expressed in terms of the surface 
pressure distributions as 
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Consequently, using (i2a) and (12b) in (9) and noting that a negative 
sign will precede the last term since Vy, 9 n < 0, 



The mathematical framework for the model of the jet-plate interaction 
j.s represented by equation 13. The physical characteristics required of 
the model can be inferred by an examination of the data already presented. 
The pertinent features of theae data are considered in the following. 

In the description of the physical nature of the jet-plate interaction, 
it will he useful to speak in general terms about a "non- destructive" 
interaction between the jet and the plate. The term non-destructive i» 
to imply that the upper portion of the jet retains its axisymmetric 
character and the influence of the plate is confined to a portion of the 
flow near the plate. The isotachs of Figures 28 to 30, the isotachs of 
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Appendix A and the photographs of the dusty jet. Figures 3b and c of [l] 
all demonstrate this effect. This non- destructive case is distinguished 
from the condition in which the approach to the plate is so abrupt that 
no portion of the jet retains its axisymmetric form. The influence of these 
two interaction conditions on the spread/curvature isotach patterns ( a - 15 
degrees) reveals that within the x/d range investigated (x/d^5) the non-de- 
structive description is applicable to all cases except a =15, h/d=l, x/d-4, 5 
and possibly a =12, h/d = 1, x/d = 5. The projected centerline of the 
jet passes through the plane of the plate at x/d = 3. 75 and 4. 75 for these 
two cases respectively. It is not possible to determine whether the non- 
destructive character is maintained for large x/d values where the projected 
centerline of the shallower angle cases would penetrate the plane of the 
plate. However, one could reasonably anticipate that for sufficiently small 
a and sufficiently large h/d, the curvature of the jet would allow the 
oblique jet to approach the character of an axisymmetric, parallel flow, 
wall jet. In this regard, the three-dimensional wall jet studies by Sforza 
and Herbst [2] are quite instructive. The upper portion of the isotach 
pattern of a rectangular wall jet does approach an axisymmetric form at 
sufficiently large downstream distances; hence, the same behavior is 
reasonable for the shallow angle, obliquely impinging jet. 

The conceptual model of the jet- plate interaction is based upon 
vorticity considerations and will be developed from the observations of 
the non-destructive interaction cases. The larger angle cases are more 
complex and additional effects must be included for their description. 

For the streamwise domain in which the jet-plate interaction is 
non-destructive the influence of the plate on the jet is manifest in two 
rather distinct regions. The phenomena occurring in these regions allow 
them to be distinguished from each other and from the approaching jet. 

In the immediate vicinity of the plate 0 Sz < 6z (where 6z/d << 1) the 
velocity distribution can apparently be characterized by isotachs which 
are nearly parallel to the plate and quite closely spaced. Such a condition 
is inferred from the isotach presentations of [l], see Appendix A, in which 
the experimental data of the y traverses indicate quite large velocities 
near the surface of the plate, and it is shown schematically in Figure 56. 
Since each lower valued isotach must lie between the measui d isotach 
location and the plane of the plate, it can be inferred t\st * v ry close 
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isotach spacing exists in this region of the flow. The region of the flow 
defined by 0 is S6z will be termed the near wall region. The second 
region lies between the near wall region and the axi symmetric portion of 
the flow. The boundary between the upper portion (of the second region) 
and the axisymmetric jet is diffuse, the two regions rather merge 
together. For convenience, this will be termed the buffer region . 

The conditions leading to the establishment of a near wall region 
are easily visualized in terms of equation 13 and the recognition that the 
isobar patterns of Figures 3 to 27 are reasonable manifestations of the 
jet-plate interaction. That is, with a pressure distribution such that 
p(x, 0, 0) forms a maximum in the p(y) distribution at a given x location, 
dp/dy will be such as to cause a flux of x-component vorticity into the 
flow. Also, the u vorticity present at the surface, as a result of the 
jet flow over the plate, will experience a du/dy reorientation effect such 
that the terms u du/dy and l/p dp/dy of (7) both cause a flux of u at 

y x 

a given x location. (The signs are -u x for y > 0 and +« x for y < 0. ) It 
is reasonable to expect that both of these effects will only influence the 
flow near the plate. The net effect is shown by the particular form of 
(13) for the spatial region defined by the limits on the integrals 
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where the u x entering the control volume at x = 0 is considered to be zero 

u du/dx << u 9 u/d u, and the signs on th«. r.h.s. terms indicate 
x y 

that negative u? x is carried out of the control volume. 

One of the contributions of equation (13) and the vortic y consider- 
ations is that the localized nature of these effects is easily accepted. 

That is, the identification of the vorticity transport effects as the 
controlling phenomenon makes the existence of a near wall region a 
reasonable result. 

Figure 56 also shows the buffer zone; this observed region is 
similarly compatible with a description of the flow in which vorticity 
transport effects govern the motion. The of the approach flow is 
essentially converted to in the near wall region; this, in turn, results 
in some production of in the near wall region as noted above. The vor- 
ticity in the axisymmetric portion remains essentially unchanged. The 
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juncture of these two regions, which occurs at 6 ^ and increasing y values 
for decreasing isotach magnitudes, results in a large and localized pro- 
duction of That is, the production term 9u/9 8 is clearly 
large at the location where the isotach makes a sharp bend to pass between 
the plate and the "undisturbed" jet. This vorticity is of opposite sign 
to that produced by the two effects operative in the nearwall region. The 
upper portion of the buffer region is not sharply defined. The isotachs 
change gradually from their circular shape in the axi symmetric region. 

This characteristic supports the implicit feature of the model in which 
the governing effects are concentrated near the surface. The upper portion 
of the buffer region is the (passive) effect caused by the vortex reorien- 
tation near the plate. 

The utilization of the vorticity transport equation has allowed 
rational explanations for the observed near wall, buffer and axisymmetric 
regions of the flow field. The same observations could be related to the 
momentum equations; however, the relationship of the acceleration, 
pressure gradient and net shear stress to the observed three regions, with 
their particular spatial extents, would have to be argued in such a manner 
that the two are simply compatible. It does not seem possible to provide 
an apriori argument for the qualitative nature of a pressure distribution 
which would create a maximum y-component velocity quite near the sur- 
face and leave the central core of the flow in an axisymmetric form. 

3.4. The Development of an Analytical Model 

In addition to its use in describing the physical nature of the jet- 
plate interaction, the model presented in 3. 3 is considered to offer an 
excellent framework for the development of a numerical computation 
scheme to describe this flow. JLighthill [3] has presented sever 1 
arguments in favor of a calculation technique based upon the governing 
equations for vorticity and not velocity. The diffusive nature of the 
vorticity equations as opposed to the action of a pressure at a distance 
allows a forward stepping soluticn to be more easily developed. It is 
anticipated that such a consideration would be most important for the 
present problem. As the jet approaches and impacts the plate, a 

description which allows the interaction to be calculated on the basis 
of the behavior at the surface and its subsequent influence on the interior 
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of the flow field is clearly superior to the alternative formulation requiring 
a construction of the three-dimensional pressure field with its consequent 
upstream effect. The development of such a numerical computation 
scheme is simply noted here; there is no immediate activity nor plan to 
pursue such a development. 


3.5. Additional Phenomena 

The principal or governing phenomena of the flow field are 
represented by the discussion in Section 3. 3; however, there are additional 
features of the flow which are considered significant. The location and 
character of the stagnation point and some quantitative aspects of the 
flux of Uy vorticity into the flow are considered in this section. 

The vorticity considerations (3. 3) tacitly assumed that the 

vorticity above the plate had a positive component, i. e. , the flow in 

the jet at (x, y, 6z) is in the streamwise direction. This is valid for the 

region downstream of the stagnation point (x g ,0,0). The stagnation point 

will occur at y=0 from the symmetry of the problem; the symmetry of the 

problem also demands that u (x.< x < x , 0, 6z)< 0 (where 6z is a small 

J ® 

z value) and x^ is the upstream boundary of the jet fluid near the plate. 

This region of reversed flow need not exist instantaneously if the stag- 
nation point moves about laterally but the mean value of u must have 
this character. The three-dimensional pattern near x g is not easily 
visualized since u > 0 can be expected at(s. < x < x g , ±6y, 6z) for 
relatively small 6y. From the isobars, it is clear that this region of 
reversed flow is relatively extensive if the stagnation point occurs at the 
maximum pressure location. The stagnation point for an obliquely 
impinging axisymmetric jet is not at the location of the maximum 
pressure; it is near the upstream edge of the interaction region defined 
by the isobar pattern. The support for this observation is given in a 
later section; its implications for the flow model are of concern here. 

The factors which cause the flow to select the location of the 
stagnation point are apparently the geometric condition that the 
stagnation streamline he located near the lower edge of the jet and the 
accommodation of the vorticity in the approach flow to that about the 
stagnation point. For example, the vorticity at the stagnation point is 
zero but its x gradient influences the orientation of the stagnation 
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streamline. Directly upstream of the stagnation point the vorticity is 
negative and of opposite sign to that of the approach flow. Directly down- 
stream, the vorticity is positive and in the same sense as that of the 
approach flow. The selection of the stagnation point is clearly a compli- 
cated physical process and it would be a difficult thing to predict. It is 
unclear what difficulties this might pose for the numerical calculation 
scheme discussed in 3. 4 in terms of formulating the forward stepping 
scheme. 

Only the x-component vorticity was considered in Section 3. 3 since 
the qualitative behavior of the flow can be described in terms of this 
variable. A quantitative description would necessitate an accurate 
assessment of as well since this component does contribute to u> x 
in the form of u^du/d y . 

The isobar patterns clearly show that the dominant vorticity flux 

into the flow is L) ; however, along the centerline only u is introduced 
x y 

through the surface. It is possible to compare the magnitude of the flux 
through the plate to the flux through a small section of 6z height at some 
x location. The analytical considerations for this are shown below. 

The flux per unit width of through an area of height 6z at any 
x location can be expressed as 

C r Sz 
\ u u dz — \ u 9 u dz 
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The flux of Uy per unit width through the plate surface from x = 0 to x = x 
is given as (from (9) and (12b)) 


x x 

^ w y dx = ' J 0 "p Tx = [p(x,0,0)-p(0,0,0)] (16) 

These two quantities can be normalized using the centerline velocity at 
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the exit plane of the jet and their ratio can be formed to demonstrate 

the relative magnitude of the vorticity flux through the surface as 

compared with that in the flow field. It should be noted that the two are 

of opposite sign for x Sx £x (max pressure) 

s 

vorticity flux through surface from 0 to x 
vorticity flux through area of height 6z at x 

2fp(0,0,0) - p(x, 0,0)1 /pn Z (o) (17) 

“u(x, 0, 6z)^/u^(o) 

Measurements of p(x, 0,0) and u (x, 0,6z) were made to evaluate 
the location of the stagnation point; the measurements and the conclusions 
regarding the stagnation point location will be presented in a later section. 
The pertinent results from this study are cited in Table 3 to evaluate the 
ratio given by (17), 


Table 3. Ratio of u vorticity flux through the plate to the 
vorticity &ux through a small height at the 
x location of the maximum pressure. Uniform 
nozzle exit condition 


a 

h/d 

x/d* 

6z/d 

{| [p(x,0,0)-p(0,0,0)} /{u(x, 0, 6z 

3 

1 

4.08 

0.0328 

14.54 / 46.2 = 0.315 

9 

1 

2. 13 

— 

27.4 / 49. 9 = 0. 35 

15 

1 

1.25 

- - - 

46.2 / 63. 6 = 0.73 


These ratios show that a significant amount of is added to the 
+u vorticity near the plate for all cases. The vorticity so added increases 
rather strongly as a increases; this is reasonable based upon the observation 
of the centerline pressure relationship to the overall and local jet curvature 
effects. The data shows that the +u y vorticity flux through 6z increases 
in magnitude even as the vorticity flux increases. This result is 
attributed to the importance of the production term due to stretching 
(u 8v/iy). Consequently the production term is quite important in the 
description of the vorticity processes occurring in the flow. Since 
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9v/9yis related to the u> x distribution, i. e. , v(x, y, 6 z) - [9 v(x, y, 0)/9 z ] 6 z = 
-u x 6 z, an analytical model would have to consider the simultaneous 
solution of the x and y components of vorticity. 


3.6. Phenomena Identified in the First Annual Report 

Several important phenomena were identified and presented in 
[l]. A summary of these phenomena is presented here for completeness. 

A control volume as shown in Figure 31, with boundaries extending 
over the exit plane of the jet and the downstream plane of constant x, the 
plane of the plate, and the sides of the flow field in the entrainment region 
will be used for this discussion. A special experiment to evaluate the 
shear force on the plate and an evaluation of the maximum probable 
pressure force acting in the negative x direction have been used to 
show that the x-component momentum flux is essentially constant for the 
oblique jet impingement problem (see [ 1 ]) . 

The impingement process has the effect of suppressing the 
entrainment of mass into the jet flow. The constant momentum flux 
combined with the decreased mass flux cause a decreased rate of energy 
dissipation. Analytically, this is described as a constant value of 




dA and a less-than normal increase 


inju 


dA which results in a 


J r 3 J 

less- than- normal decrease in \ u dA, where the integration is over a 
plane normal to the axis of the jet. Mechanistically, it is inferred that 
the decreasing energy flux is a measure of the work done by the jet fluid 
to accelerate the ambient fluid as it is entrained into the jet flow. Since 
the plate inhibits this entrainment process, the relatively high velocity fluid 
in the jet is protected from the decelerating effects of the entrainment 
process and the only energy decrement incurred in this region is that as- 
sociated with the boundary layer in the near wall region (0 Sz £ 6 z). The mass 
and energy flux values are shown in Appendix C. The data for the 
uniform and fully developed nozzle conditions are shown separately 
and the data are segregated with respect to how well the measured 
x-component momentum flux agrees with the quantity J(o)cos 9. 

The data for which the measured momentum flux is within three percent 
of that at the exit plane is considered to be quantitatively valid. That which 
differs by more than three percent from the constant value is considered 
to be qualitatively instructive in that it shows the appropriate trends for 
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tin* mass and energy flux values. 

An interesting difference between the uniform and fully developed 
exit plane conditions is the greater rate of entrainment of ambient fluid 
for the fully developed condition. It is inferred that this greater entrain- 
ment rate is related to the large scale motions in the initial region of 
this flow. The mass rates of flow for the two nozzles of the present 
study and for several other studies which have been reported are also 
presented in Appendix C. A comprehensive study of the initial condition 
effects on the near field of an axisymmetric jet is presently being 
conducted as a part of the continuing activity under this grant. 

The jet- plate interaction can be considered complete in some sense 

when the total z-component momentum of the initial jet is destroyed. 

Therefore, the degree of compLteness of the interaction is shown in 

terms of the ratio of the net force on the plate to the initial z-component 

,x p® 

momentum flux. Curves of \ \ p dy dx/j (o) are presented in [ 1 ] . 

^0 J -® z 

Alternatively (and less graphically) the z-component momentum flux, 
normalized by the initial z-component momentum flux, provides a measure 
of the degree of completeness of the interaction. The tabular 
representation of these values taken directly from [1] is presented in 
Appendix D. 

4. THE CENTERLINE PRESSURE DISTRIBUTIONS 

From a mechanistic viewpoint, the surface pressure is the agent 
which causes the spreading and jet curvature imposed upon the jet by the 
no- penetration condition through the plane of the plate. On this basis, 
the only apriori constraint on the surface pressures is that the pressure 
integral over the plate surface is equal to the z-component momentum flux 
from the jet. One aspect of the total pressure distribution is the pres- 
sure along the line (x, 0,0); that is, the centerline pressure distribution. 

In Section 3. 2 the centerline pressures were compared with the curvature 
of the jet to demonstrate the local, in addition to the global, effects of 
the curvature. The essential point is that the centerline pressure distri- 
bution is governed by the geometric and, to some degree, the initial flow 
structure conditions of the jet. In this sense it is an effect, the no- 
penetration condition is the cause. 
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A different condition exists for the streamwise vorticity. The 
centerline pressure can be interpreted as the cause of the amount of x- 
component vorticity added to the flow by the "flux through the plate surface. " 
Equation 18 provides the basis for this statement; specifically. 


u x introduced through the 

differential area dxdy via 
the pressure distribution 


8 u 

VJ t~ d y dx 


■ p (18) 


The flux through a given area of the plate cannot be directly related to the 
flux through a specific y-z area at some x location. The reason for this is 
that the velocity field rear the plate will convect and the velocity gradients 
will amplify or reorient vorticity which enters the flow field through this 
mechanism. However, it is possible to identify the total vorticity which 
enters the flow and to attach some meaning to this measure of the u> x 
vorticity flux. Consequently, (18) is integrated over the region from 
y=0 to + (or -) oo and from 0 to x. The resulting value is that total 
<d x introduced into the flow via this mechanism; specifically 


id introduced via the pressure 
x 

distribution 



— dy dx 

p 8y 1 



-■ [p(x,«iO) - p(x, 0,0)] dx 

P( X *°>.QL dx 
P 


(19) 


Equation 19 demonstrates the importance of the centerline pressure 
distribtuion to the flux of These integrated values are presented in 

a composite form in Figure 57. Only the fully developed cases are 
shown for a < 15 degrees, the uniform ones are essentially the same. 

The strong dependence of the vorticity flux upon a , and h/d for 
a given a , is evident. The strong differences shown for the various 
a and h/d cases suggest that the flows near the plate surface are 
strongly dependent upon these two parameters. However, the data of 
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Table 2 show that the width of the flow field and the pressures in the 
neighborhood of the 0. 1 free isotach intersection are weakly dependent 
upon these two parameters for the non- destructive interaction cases. 
Therefore, it is inferred that the differences in the flow structure 
associated with the flux of u> x are confined to the central region of the 
flow field. 

5. THE STAGNATION POINT 
5.1. General Considerations 

An important characteristic of the oblique jet impingement which 
must be established for a complete description of the flow field is the 
location of, and local behavior around, the stagnation point. The 
identification of the stagnation point characteristic is considered to be of 
particular importance for those characteristics which depend upon the 
details of the flow field near the surface such as the heat transfer, 
acoustic noise generation, and the production of streamwise vorticity. 

The stagnation streamline is defined as the locus of points 
everywhere tangent to the velocity field with a terminus on the impact 
plate. The stagnation point is the location of the terminal point. The 
possibility of a finite length stagnation line will be suggested by the 
experimental data and it cannot be excluded by analytical considerations. 
The possibility of a finite length stagnation line is also suggested by the 
following considerations. An infinitely long stagnation line exists in a 
two-dimensional flow and a stagnation point or circular stagnation 
line exists in an axisymmetric flow. Similarly, an open conduit of any 
shape placed with its axis parallel to a uniform stream will result in a 
closed, stagnation line around its upstream circumference. 

The three-dimensional flow may be qualitatively different from the 
axisymmetric or two-dimensional condition. For example, consider the 
open conduit placed in a uniform stream and consider a sawcut to be 
made through the surface parallel to the axis. If the only possible 
stagnation states are a point or a closed line curve, then a single 
stagnation point would have to exist opposite the sawcut. Since the flow 
will still pass inside and outside the contour defined by the previously 
closed body, the symmetry would demand that an azimuthal flow exist 
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from the stagnation point, around the upstream surface and exhaust down 
the gap of the sawcut. The alternative is to conceptualize a three- 
dimensional stagnation line extending part way around the upstream surface. 
This aside from the oblique jet impingement problem is to demonstrate, 
by analogy, the possibility of a finite length stagnation line. The oblique 
jet impingement is a three-dimensional flow. The possibility of a finite 
length, non- closed, stagnation line would seem to exist. 

Another characteristic of the stagnation point is its relationship 
to the location of the maximum pressure point. For a vertical 
(a = 90 degrees) jet impingement condition, it is clear that a stagnation 
point exists at the terminus of the center streamline and at the location 
of the maximum pressure. The centerline pressure data shown in the 
previously discussed figures indicate that the maximum plate pressure is 
not at the projected location of the center streamline; these figures also 
suggest that the maximum pressure may not be located at the stagnation 
point. The analytical considerations which can be made for the stagnation 
point are presented in the next subsection. An experimental technique 
to evaluate the stagnation condition and the results of this evaluation 
for three cases are presented in this section. 

5.2. Analytical Considerations* 

The equations of motion (Navier-Stokes and continuity) can be 
evaluated in the neighborhood of the stagnation point in order to 
establish certain characteristics of the stagnation streamline. These 
analytical considerations are presented in this subsection (5.2). 

5.2.1. The orientation of the stagnation streamline 

Let (6x, 0, 5z) be the coordinates of a position on the stagnation 
streamline. The stagnation point will be designated (0,0,0) for the 
purposes of this analysis and the symbol 5s will designate an incremental 

*The remainder of Section 5 is adapted from the second semi annual 
report [4], A portion of the analytical considerations was incomplete 
although the results were not affected and the experimental data have been 
more completely analyzed. For these reasons, and because of the limited 
distribution of the semi annual report, these subsections are repeated 
herein. 


25 



2 2 1 /2 

length along the streamline where 6s = [6x + 6z ] 1 and 

6x/d~ 6z/d << 1 is understood. The slope of the stagnation streamline 

at the stagnation point is given as 

5 7 * 

slope of stagnation = lim 7- (20) 

streamline at (0,0,0) 6s-* 0 x 


and this slope can be expressed in terms of the velocity field as 


lim 
6s- 0 


6z 

6x 


lim 
6s— 0 


w(6x, 0, 6z) 
u(6x, 0, 6z) 


( 21 ) 


Assuming that the u and w components are analytic functions , the limit of 
the velocity ratio can be expressed in terms of the Taylor expansion of 
the velocities. Only the lowest order, non-zero terms of the expansion 
need be retained in anticipation A taking the limit. All partial 
derivatives are understood to be evaluated at (0,0,0); 9u/9x = 9 w/ 9x =0 
because of the no slip condition and 9w/9z = 0 by continuity. 


wj6x , ,0 J _6zj 

6x— 0 u(6x '°' 6z) 


6z 

2 


9 2 w 

a 2 
9 z 


6z 


9 u 
9 z 


6 z 

2 


9 2 u 

9z 2 


+ 6x6 z 


a 2 

9 u_ 

& xd z 


( 22 ) 


= 0 for 9 u/9 z / 0 
>, < ,= 0 for 9 u/9 z = 0 


(22a) 

(22b) 


Condition (22a) shows that if the vorticity at a point is non- zero then the 
streamline at that point (in the direction perpendicular to the vorticity) 
lies in the plane of the plate and hence cannot be a stagnation streamline. 
Equation 22 can be operated upon to provide a more easily interpreted 
form when (22b) is applicable. Dividing by 6x6z yields 


lim 6z = 
6 s- 0 6x 


i. 5 z 9 2 w 

2 6x a 2 ‘ ' • 

o z 

J. 2 JL 9 2 u . 9 2 u 

2 6x . 2 9 x9 z 

9 x 


(23) 
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The second derivative of u with respect to x and z can be expressed in 
terms of w from the continuity equation as 



8 2 u 8 8 v 8 2 w 

” 8x8 z 8z8y 2 

OZ 

and since d v/dzdy = 0 by symmetry, 


a 2 A 2 

d u_ _ _ 8 w 

dxdz a 2 

O 2 


Consequently (23) can be written as 


lim 
6 s -0 


6 z 
6x 


1 

2 


,1 

2 


6 z 8_jw 
6X 8 zF 
5z 8 2 u 
6x 8x 2 



(24) 


(25) 


Equation 25 can be examined for the orientation of the stagnation 
streamline for the condition that (1) the stagnation point occurs at the 
maximum pressure point and (2) that the stagnation point occurs at a 
location where the streamwise pressure gradient is positive. The 
8 u/8 z c derivative can be related to the pressure gradient by the 
x-component momentum equation evaluated in the plane of the surface; 
specifically, 


8 p 8 Z u 

^ = v 77 


for z = 0 


2/2 

Consequently, for the condition 8“u/3z ]q 
25 becomes 


(26) 

8p/8x ]q = 0, equation 


lim 

6 x -*0 

6 z -*0 


6z 

6x 


i. §5. 

2 6x 


for 


i£_ 

8 x 


0 


0 


(27a) 


which is only possible if lim 5z/6x = 0 or ®. A zero value cannot 


27 



represent a stagnation point; hence, the stagnation streamline intersects 
the plate at an angle of ir/2 when the stagnation point occurs at the 
maximum pressure. This would occur for the a = it / 2 condition. 

For dp/dx / 0, equation 25 can be expressed as 


lim 

6s-0 


6z 

6x 


6z 

6x 


9 Z u 
8 z 2 
8 2 w 


8 z 



2 


(27b) 


and (27b) shows that (note that the denominator must be equal to unity) 


*2 
o w 

6 z _ 8 z 2 

6x a 2 

9 u 


(28) 


Using (24) and noting that 




(29) 


The streamline orientation at the stagnation point is 


lim 
6 s— 0 



-v 


1 

P 



8 x 

li> 

8 x 


(30) 


That is, the orientation of the stagnation streamline earn be expressed in 
terms of the pressure and surface vorticity gradients at the stagnation 
point. 

5.2.2. Consequences of stagnation point at or displaced from the 
maximum surface pressure 

The relationship between the vorticity which is introduced into 
the flow field by the viscous shear at the wall and the static pressure on 
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the plate has been presented in Section 3. 3, viz. , 

(12b) 

(12a) 

z .-0 

Since the solenoidal condition on vorticity demands that it appear in closed 
loops (i. e. , V • u, = 0), it is of interest to note the relationship between 
the created vorticity and the stagnation point considerations. 

Consider the case of vertical jet impingement. In this condition, 
the stagnation point is centered in the axisymmetric flow field and the 
stagnation point is clearly at the maximum pressure. That is 
lim 6z -*■ 0 6x -*■ 0 w/u - -<^ and 9 p/9 x = 0 = 9 p/9 y. A fluid element 
immediately adjacent to the stagnation streamline will possess vorticity 
in the sense indicated by Figure 58 . The vorticity of such an element must 
change sign along its trajectory from above the plate to a location in the 
boundary layer flow along the plate. This vorticity sign change is accomplished 
by the viscous diffusion term of the vorticity transport equation, viz. , 
from (8) 

=Zj° V & + i H + u * Vd - 6 • V (b (8) 

Dt 

(1) (2) (i) (4) (5) 

(1) convective transport and local time rate of change of vorticity, 

(2) ’'production" of vorticity by vortex filament stretching and re- 
orientation, 

(3) diffusion of vorticity, 

(4) production of vorticity from the fluctuating vorticity- velocity 
interaction, 

(5) vorticity transport by the turbulent motion. 

As shown by Figure 58, the flux of vorticity at z=0 is symmetric 
about the stagnation point, that is, the new vorticity introduced into the 
flow is created in symmetric loops about the stagnation point. Also, the 
fluid which originally comprised a vortex loop will remain in the same 
loop (on the time average) as the sign of the vorticity is changed along the 
trajectory. This rather straight forward condition provides a useful 
reference for th.e considerations of oblique jet impingement. 


i i£_| - jl 

M ax dz 


2 y 


r. = a 


and ± |E 
U oy 


z=0 


7T- W 
9 z x 
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For an oblique impingement, the center streamline will not, in 
general, be the stagnation streamline; hence the vorticity structure near 
the stagnation point will be quite different from the normal impingement 
case. This is shown schematically in Figure 59. For example, 
the fluid in the original loop 1-la is not in a common vortex loop 
following the impingement; that is, the vorticity of the fluid at 1 changes 
sign, the vorticity of the fluid at la does not. Consequently, the vortex 
and material lines are different for the oblique jet impingement case. The 
£ - 0 surface will be non- symmetric and strongly three-dimensional for 
this flow. It is not obvious what shape this surface will assume away 
from the centerplane. Smoke flow visualization studies have shown tnat 
the entrained flow stagnates and then separates from the surface as 
shown in the sketch. 

It is not clear what geometric configuration is assumed by the 
closed vortex loop involving the fluid at 1 when it has progressed to a 
position just upstream of the stagnation point. It seems likely that these 
loops which are created following the stagnation process on the adjacent 
streamline are somewhat symmetric about the stagnation point. In this 
regard, the production of the new vorticity would be similar to the 
processes occurring in the normal impingement case. 

5.2. 3. .Relationship between the stagnation and surface pressure gradients 

near the stagnation streamline 

An analysis of the surface and stagnation pressure gradients will 
be necessary for the proper interpretation of the experimental 
measurements. Specifically, we will show that the surface pressure in the 
backflow region near the stagnation point is less than the stagnation 
pressure of the streamline in the neighborhood of and located upstream of 
the stagnation streamline. Figure 5V shows the terms for the analysis. 
Point A lies along the stagnation streamline point B is the intersection 
of the normal at A and the second streamline of interest and point C is 
at the intersection of the second streamline and the contour which lies 
along the normal trajectory of the streamlines from x g - $ x to the 
streamline passing through B. 
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The pressure at (x - 6x, 0, 0) can be written as 
s 


p(x g - fix, 0,0) = p(c) + ^ dn 


( 31 ) 


where L is the distance from C to (x - fix, 0, 0). The stagnation pressure 

S 

and the velocity at C can be used to express p(c) as 


p(c) = pJ c ) - ^ 


(32) 


and the integral along n can be written as 


$0* d " 


c 1 pQ 2 dn 

J 0 B/L) L 


(33) 


The maximum modulus theorem can be used to express this integral 
value as 


f l ^a! *> < 

Jorn/L) L 


pQ 2 ) 


max 


R/L) . 
min 


and since Q 


max 


>o(RA) 

Q(c) equation (26) becomes 


p(x s -6x,0,0) < p s (c) - 


min 


(34) 


(35a) 


and 


p(x -fix, 0,0) < p (c) if R/L > 2. 
s s 


(35b) 


The dimension I, is related to the displacement of B with respect to A 
and consequently can be made as small as desired. Hence, the condition 
R > 2L is assured, since R is finite as L-* 0. (Note, R -*■ ® as L *♦ 0. ) 

It is now necessary to show that this result is general; that is, 
that it applies for all small 6z. For sufficiently small L, the stagnation 
pressure at C (or B) and the pressure at (x g -6x, 0,0) can be related 
to the values at A and (x ,0,0) respectively using a first order Taylor 
series. Hence the inequality of (28) can be used to express an inequality 
between these pressure gradients. Specifically, assuming the stagnation 
pressure changes to be negligible over these small displacements. 
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w 

(37) 

(38) 

(39) 


geometrically, it can be shown as in Figure 59 that 6x > 6n for small 
6x. Consequently 3p/3x > 9p s /9n. 

The important result of these considerations is that the surface 
pressure is always less than the equivalent stagnation pressure except 
at the stagnation point. This is a necessary result for the proper 
interpretation of the stagnation point. 

5. 3. Experimental Investigation of the Stagnation Phenomena 
5. 3. 1 . Technique and interpretation 

Stagnation phenomena for the oblique jet are difficult to investigate 
because the phenomena are basically unsteady and hence a suitable 
averaging process is required. Various schemes to determine the time- 
mean stagnation streamline by use of its singular characteristic as the only 
streamline to intersect the plate were unsuccessfully investigated. Holographic 
interferometry with a dense gas tracer and a heated jet with a platinum 
resistance temperature sensor (high frequency response = 3khz) were 
used to attempt this discrimination. The results from the attempts were 
not decisive. The common and often successful technique of attempting 
to locate the position of zero shear stress was not attempted because the 
expected stress magnitudes for these small angles would be quite small 
and hence difficult to discriminate and because this flow is strongly 
three-dimensional which would necessitate a measurement over a small 
length and width ctmpounding the sensitivity requirements. Also, a suitable 
averaging would have to be inherent in the transducer or in the subsequent 
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signal processing in order that the time variation would be properly accounted 
for. 

An experimental technique involving simultaneous hot-wire and surface 

static pressure distributions was developed and successfully employed. 

For this technique, the hot-wire probe was positioned close to the plate 

(zS30. 15 cm, z/d = C. 04) and the velocity magnitude was recorded as a 

function of the distance from the jet v w (x). A pressure survey for the 

same conditions allowed the calculation of a velocity as if the surface 

pressure was the stagnation pressure for that streamline, this value was 

designated v . A linear shift 6x of the hot-wire velocity data to the z-value 

of the plate, along the line defined by the wire position z and the jet angle 

a plus the divergence angle (3 of the 0. 1 isotach, was then used to provide 

a comparison between v and v . The equation for this shift is 
_ j P w 

Ax = Az tan (a + (3). The stagnation point is considered to he at the 
intersection of these two curves. A sketch of this technique is shown in 
Figure 60. 

5. 3. 2. Experimental results 

Three cases were selected for comprehensive evaluation of the 
stagnation point location, h/d = 1 and a = 3, 6, 9 degrees for the"uniform 


nozzle" configuration. The results of the special traverses to determine 

v and v are shown in Figures 61 to 63. . Li general, two v traverses 
w p 6 w 

were made at different z/d values. The results from the twj traverses were in 

reasonable agreement for the three different angles; the meet complete 

set (i. e. , z/d value) was used for these results. The data themselves 

are presented in Table 4. The results for these three conditions show 

that the stagnation point does not occur at the location of the maximum 

surface pressure. The v and v data for a = 3 degrees, h/d = 1, 

■\* P 

(see Figure 61), did not indicate that a stagnation point existed. However, 

when the velocity data from the hot-wire probe were shifted an additional 

distance of 0. 57d a stagnation region is evidenced. The requirad shift is 

interpreted as an upward deflection of the jet by the increased static 

pressure on the plate. The location of the jet's center, based upon 

momentum flux considerations (z of the first annual report), dees not 

m 

show a shift above the geometric trajectory for the condition a = 3 degrees, 
h/d = 1, uniform; such shifts are indicated for a = 3 degrees and h/d = 

0. 75 and 0. 5. However, it is reasonable that the low velocity fluid along 
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the stagnation streamline could be affected by the surface pressure for the 
h/d = 1 case. 

The location of the stagnation point is quite difficult to infer from 
these results. Apparently a more accurate description is to identify a 
stagnation region in the neighborhood of x/d = 3. 8. The flow field is 
considered to possess a single stagnation streamline at an instant; however, 
the inherent unsteadiness of the flow could well produce an apparent stagna- 
tion region even though the time mean flow field will also possess, by assump- 
tion, a unique stagnation streamline. If viscous effects reduced the stagna- 
tion pressure between the velocity measurement location and the surface, 
then this additional shift would not be as large as is indicated. A shift of 
Ax/d =1.5 for the hot-wire velocity data ( v w ) would cause the stagnation 
and maximum surface pressure locations to be coincident. However, 
this would result in a condition in which the local surface pressure would 
exceed the equivalent stagnation pressure; this condition is shown to represent 
a physically impossible situation in equation 39. It should be noted that 
a linear shift of the data over a net distance of Ax/d = 1. 20 is too great 
to allow detailed comparison between the v^ and results. 

The data for the case a = 9 degrees, h/d = 1, u, see Figure 62, 
appear to be quite straightforward in terms of locating the stagnation 
point, x/d = 1.9, and demonstrating chat the stagnation point does not 
occur at the maximum pressure location. The absence of a vertical 
deflection is compatible with the z^ and z m results for this case (see Figure 
45d, First Annual Report) and represents a reasonable trend as compared 
with the three degree case since the jet approaches the plate with a less 
glancing trajectory. A slight negative shift is suggested (but not shown) to 
provide a tangential condition between the v^ and v^ curves. 

The measured v and the measured- shifted v distributions for 
p w 

a =15 degrees are shown in Figure 6 3. These results suggest that the 

jet trajectory is rteeper than 19 degrees in the region between the plane 

of the measured velocities (z/d = 0.0308) and the surface; hence, the 

indicated shift is too great. This observation is based upon equation 39 

in which it was shown that the surface pressure should not exceed the 

equivalent stagnation pressure. The solid line on Figure 6 3 presents the 

estimated upstream shift appropriate for this case. The resulting 

tangent poir.t between v and v i«» in the neighborhood of x/d = 1. 3. 

p w 
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5. 3. 3. Summary and analysis of the stagnation point investigation 

The determination of the stagnation point by the special traverses 
discussed above allow the relationship between the general velocity field 
and the stagnation point to be examined. The approximate locations of the 
stagnation points for the three cases are 

uniform, h/d = 1, a =3 x/d = 3.8 

9 1.9 

15 1. 3 

A comparison of these results to the isotach plots of Figures 28 to 30 
shows that the stagnation point is in the vicinity of the intersection 
of the jet edge (u/u(0) = 0. 1 isobar) with the plate. 

There are two physical effects which could cause the inferred 
stagnation point to be in error; they are self-compensating and their 
magnitudes are difficult to estimate although they are certainly too small 
to affect the qualitative nature of the above conclusions. Shear effects 
will degrade the magnitude of the stagnation pressure between the measure- 
ment location and the surface of the plate. To the extent that this occurs, 
the measured velocity represents an excessive stagnation velocity and the 
v^ magnitudes would be corrected by subtracting a prescribed 6v w value. 
The physical proximity of the velocity measurement and the surface 
pressure measurement would cause this to be a relatively small error. 

If the velocity were measured at a z location where p> P arn ^ e then 
the stagnation pressure is larger than that indicated by the pv^/z value. 
Hence a correction of adding a separately prescribed 6v^ would be required 
to relate the hot-wire and surface pressure data. It is difficult to estimate 
reasonable values for for the static pressure correction; however, it 
can reasonably be assumed that this correction increases as the wall 
pressure increases, that is, d(6v w )/dx > 0. Consequently, the alignment 
of the v^ and v^ curves would always occur at a point farther upstream, 
allowing for this correction, and this reinforces the result that the 
maximum pressure and the stagnation point are physically displaced. 

6. SUMMARY 

The kinematic, no penetration constraint requires that the fluid in 
the jet spread laterally or deflect vertically upward to provide the required 
cross sectional area for the jet mass flux. 
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The selection between these two responses is based upon the 
dynamics of the jet-plate interaction as characterized by the momentum 
equation. 

Analytical and experimental considerations are used to demonstrate 
that for sufficiently small angles, the primary response of the jet is an 
upward curvature; however, a complex interaction between curvature and 
spreading is suggested for angles larger than (nominally) 12 degrees as 
would be encountered in the interaction between the jet and the deflecting 
flaps of a STOL configuration. 

The description of the jet-plate interaction in terms of vorticity 
is well motivated by both analytical and experimental considerations. 

Such a description provides a rational framework to model the three- 
dimensional character of the flow in which the upper portion retains its 
axisymmetric structure and the distortion of the jet occurs primarily near 
the surface. Streamwise vorticity is introduced into the flow field by the 
lateral surface pressure gradient and the reorientation of the azimuthal 
vorticity of the approaching jet. 

Similarly, there is a surface flux and a production by stretching 
of the lateral vorticity. Such phenomena account for the localized jet 
distortion. The streamwise vorticity could be used to provide enhanced 
boundary layer control if properly accounted for in the flap gap design 
for an externally blown flap airfoil. 

The physical displacement between the stagnation point and the 
maximum pressure point is examined analytically and demonstrated from 
the experimental data for incidence angles of 15 degrees and less. The 
possibility of a finite length stagnation line is demonstrated for a non- 
analytic velocity distribution and the relationship between the incidence 
angle at impact and the surface vorticity and pressure gradients is 
established assuming an analytic velocity field. 
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APPENDIX A. 


Selected Isotachs from the First Annual Report 
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3 2 1 0 -1 -2 -3 

a) x/d * 0 


Figure 14. Isotach contours a = 3, h/d = 0. 75, fully developed. 
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c) x/d = 2 
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d) x/d = 3 
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e) x/d = 4 
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i) x/d = 5 
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3 2 1 0 -1 -2 -3 

a) x/d = 0 


Figure 18. Isotach contours a = 6, h/d = 0. 75, uniform. 
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3 2 1 0 -1 -2 

b) x/d = 1 
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c) x/d = 2 
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3 2 1 0 -1 -2 -3 

a) x/d * 0 

Figure 21. Isotach contours a = 9, h/d = 1, fully developed. 
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b) x/d = 1 
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3 2 1 0 -1 -2 -3 

c) x/d = 2 
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3 2 10 - 1 - 2-3 

a) x/d * 0 

Figure 22. Isotach contours a = 12, h/d = 1, fully developed. 
















3 2 1 0 -1 -2 -3 

b) x/d = 1 
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d) x/d = 3 
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a) x/d = 0 

Figure 24. Isotach contours a = 15, h/d = 1, fully developed. 
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d) x/d = 3 


64 











Symbol 

o 

□ 

□ 

□ 

a 

Q 

VVA 

□ 

□ 

A 

o 

0 

u/u(0) 

m 

3H 

0.9 

0.8 

SB 

QQ 

m 


0. 3 

m 

0. 1 

0. 05 


h(x) = h - x tan a 

s = s location of the 

_ _ 

_ isotach-circle center 

i 

o 

- o „ 

o 

o 

o 

o 

o a 

o' 

1 

□ 

> 

□ 

o 

7 U u 

0 0 n cP v7 A - 

□ V A 

- ^ °° 

• 

A 

◄ 

□> 

o 

v*. - 

1 J n 1 — ▲ 1 • — 

1 A 1 L_ 


3 2 10-1-2-3 

e) x/d = 4 


65 













Symbol 

□ 

□ 

a 

D 

□ 

u 

VA 

□ 

□ 

□ 

D 

a 

u/u(0) 

0. 99 

m 

SI 

0.8 

QQ 

m 


SI 

SI 

Si 

0. 1 

0. 05 


h(x) = h - x tan a 


* c = * location of the 

iaotach-circle center 
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Figure 25. Isotach contours a = 15, h/d = 1. 5, fully developed. 
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3 2 1 0 -1 -2 -3 

b) x/d = 1 
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APPENDIX E. Corrections to the Center of Momentum z Analysis, 
Section 4. 5 . 2 of [l ] m 

Following equation 30, the expression v = - n Q results in 

c c 

(v e • k) (v e • n) = + n • k (A. 1) 

and hence equation (31) rhould be (the c following the equation number 
denotes "corrected") 


F (x) , n 

t * pdA 


J 

jftj < x > + sinQ 


2 J(0) 


1 


pQ 


A 

k 


A 

n 


dA 


(31c) 


The effect of the low pressure in the entrainment region was 
omitted from the moment- of- momentum equation (32). Including this term 
results in the following expression for the negative y-component equation 
evaluated with r measured from th*- '*er of the exit nozzle 



x dA 


(h - 


z ) J 
m x 


+ xJ 


\ p(rxv ) 
J A 


— A I A 

• j(v - 1 , oA - \p(r x n) 

J A 

e 


A 

« j dA 


(32c) 


The remark that "the entrainment integral is zero except for the 
plane x = 0" is incorrect. Consequently, this term mu6t be included in 
equation 33. The integrand of the entrainment integral may be written 

<L j 


P I(r x v e ) 


J] l(v 


A , 

n ) J 


p{-(z n • i - x n • k)Q e J [- Q e J 

A A A A t 

p [z n • i - x n ® k] Q g (A. 2) 


and combined with the integrand of the pressure integral over the 
entrainment area as shown by the following. The integrand of the pressure 
integral is 

A x A A A A A 

P e <r x n) o j = P e lz n • i • x n • k] 
or, using (29) (i.e., (p ) = - pQ^/2) 

eg e 
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- A a °e r A a a \ 

p(r x n) « j -p — : — |z n • 1 - x n • k| 

Ct 


and combining this with the integrand of the other entrainment integral of 
(32c) results in 


l C 0 2 r A A A (') 

- | J pQ^ [z n • 1 - x n • kj dA 


(A. 3) 


The corrected form for equation (33) can now be developed 


m 


, J 

h x z 

d d J( 0 )cosq 


1 f X 1 C A A A a 1aa 

iioicos; j/d 44 • 'm>cosa}p Q J z r •*-* n * k i dA 

P e 


h x p x x 

" d d J(0) cosa " d tan a 2J(0)cosa 


/d C ? A A 

J A " Q e K -" dA 


J(0)c 


- h + 2£ £ 

1 1 IT /A \ 


cosa d P d ^ 2J(0)cos a ^ 
P 


? A A A A 
pQ g (z n « i - x n « k) dA 


d d J(0)cosa d 


- ^ tana 


i— r $ 


J(0)cosa d 

P 


p dA 


1 


2J(0)cosa ^ 5* P °e * d “ * * " d " •' k)dA + d i/ k * " dA S 
L e e 


(33c) 

Reasonable approximations can be made to this moment- of- momentum 

equation. Specifically, in the region where the jet is curving because of 

A A 

the plate presence, the term z n « i can be considered small with respect 
A A 

to the term x n • k. Consequently, the bracketed term in the above equation 
can be approximated as 

1 If rt Z,z A * X A , . . x f ^2 A . A , .1 

_____ pQ( — n«i--rn»k) dA + -r\pQ k°n dA r 

2J(0)cosa |~A e d d d J A e j 

rj P Q e f " • * dA + 2 Q“ n . k dA \ 

l A p J 


2J(0)cosa 
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(A. 7) 


1 

l J(0)cosa 



A 

k dA 


/' further simplification is possible on the basis of the numerical 
magnitude of this term with respect to * m /d. The jet will be deflected 
upward in such a manner that the z m /d values will be of the order 
0.25 or greater. (See the calculations of the First Annual Report [l] 
which, although in error, are of the correct order of magnitude. ) The ratio 
of the entrainment velocity at the edge of the jet to the jet exit velocity 
can be evaluated from the mass flux evaluations of the First Annual 
Report [ 1 J. In order to estimate the older of magnitude of the 
entrainment integral term we use the result for the value of 
d[M/(K f)^^pu(o) A Q ]/dx from Figure 34 of £ e earlier report. The 
values are 

[ M/(k \)*'^pu(o) A q ] = 0.171 fully developed nozzle (A. 4) 
pu(o) A q ] = 0.134 uniform nozzle (A. 5) 

where k and \ are coefficients to correct the momentum flux and to 
allow u(o) to be used as a normalizing velocity respectively, see [l]. 

The value of k can be considered to be unity. The velocity magnitude 

O 

Q e can be related to the M values by the relationship 

p 2ir R Q = -j— & \A.6) 

r e dx 


where R is the radius at which Q is evaluated. Since the ratio 

e 

2J(0)Lr I P °e a , k jA is desired it is only important to accurately 

A A 

model Q g where the product n • k is appreciable. This is fortunate since 
the isotachs of the First Annual Report indicate that the ;pper portion of 


the jet is unaffected t>y the jet-plate interaction except for a sufficiently 
large combination of a and x/d and even then the upper portion of the 
jet tends to retain a circular appearance. Also, since the control volume 


is a semi circular arc above the plate and since the entrainment in the 

A A 

bottom of the jet enters the lateral sides where n • k = 0 it is only 


necessary to accurately model the upper portion. The magnitude of 

• . 

can then be evaluated from (A. 6) as (note that the linear dM /dx 
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implies dQ^/dx = 0) 




Recognizing that (d/2R) < 1, the magnitude of the term in equation 

(A. 7) can be estimated as 


1 

2J(0)cosa 



= (V 1 ) khi) <°- 7 > (I) (f ) 2 

s (x/d) 2 x 10" 3 (A. 9) 

A A 

where (0. 7^ is considered to be the average value of n * k and the 

. 3 

combination of the rumerical values is rounded to 10 

These considerations show that the entrainment effects are 
small for the x/d v 0 s x/ d ^ 5) domain of interest. Hence, the integral 
representing the entrainment effects can safely be neglected for the 
purposes of computing z (x) and the quantities which depend upon it. 
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APPENDIX C 


Mass, Momentum, and Energy Flux Values 
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Figure 30. Normalized mass and energy flux values, uniform, 0 S | k - l| S 0. 03 
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Figure 31. Normalized mass and energy flux values, uniform, | k - ij >0. 03. 
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Comments The data of Hill 1 19], obtained with 
a porous cylinder, indicate that the slope of the 
mas'j flux curve is not constant in this region. 
The relative accuracy of the porous cylinder 
versus the integration of the velocity readings 



Figure 1. Normalized mass flux values as evaluated from references 
[12], { I 4 j, and [2l] (Data for Crow and Champagne [29] not 
not shown, a = 0. 1 36 for 0 £ x/d 5 2; a j = 0. 292 for x/d > 6, ) 
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APPENDIX D. 


Normalized J Values to Demonstrate the Effective Extent 
z 

of the Interaction Region 
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- 0.99999 


0.99999 


Z -COMPONENT MOMENTUM/ SI N< ALPHA) 

15.000 

12.000 - 1.00000 

ALPHA 9.000 - 0.99999 

6.000 

3.000 

x/i = 1 0.000 


x/d = 2 


x/d = 3 


x/d = 4 


x/d = 5 



0.50000 

0.75000 1.00000 

1.50000 

2.00000 




H/D 


Z-COMPONENT MOMENTUM/SINIALPHAI 



15.000 


- 0.80419 

- 

1.00529 

12.000 



- 1.00951 


ALPHA 9.000 


- 0.99999 

- 1.00192 


6.000 



- 0.99999 


3.000 





0.000 






0.50000 

0.75000 1.00000 

1.50000 

2.00000 




H/D 


Z-COMPONENT MOMENTUM/SINIALPHAI 



15.000 


- 0.39684 

- 

1.01222 

12.000 



- 1.00922 


ALPHA 9.000 


- 0.86540 

- 1.00471 


6.000 


- 0.99999 - 0.99999 

- 0.99938 


3.000 


- 1.00000 - 1.00000 



0.000 


0.00000 0.00000 




0.50000 

0.75000 1.00000 

1.50000 

2.00000 




H/D 


Z-COMPONENT MOM ENTUM/SINI ALPHA) 



15.000 


- 0.10517 



12.000 



- 0.92929 


ALPHA 9.000 


- 0.61194 

- 1.00707 


6.000 


- 0.60765 - 0. 9245 

- 0.99814 


3.000 


- 0.74821 - 0.96427 



0.000 


0.00000 0.00000 




0.50000 

0.75000 ..00000 

1.50000 

2.00000 




H/D 


Z-COMPONENT MOMENTUM/SINIALPHAI 



15.000 





12.000 



- 0.68396 


ALPHA 9.000 



- 0.94577 


6.000 


- 0.38323 - 1.23213 

- 1.00061 


3.000 


- 0.50039 - 0.87497 



0.000 


0.00000 0.00000 




0.50000 

0.75000 1.00000 

1.50000 

2.00000 


M/0 


C 2. 


Z-COMPONENT MOMENTUM/ SI Nl ALPHA) 

15.000 

12.000 

ALPHA 9.000 

6.000 

. 3.000 

x/d = 6 0.000 0.00000 


1.41309 - 0.97193 
0.73271 


0.50000 0.75000 1.00000 1.50000 

H/P 


2.00000 


Table l.a. Normalized J values to show the effective development 

length for the uniform nozzle exit condition (J z /J(0) sin a). 
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x/d 
x/d - 
x/d, 
x/d = 

x/d = 
x/d = 

Table 1.1 


Z -COMPONENT MOMENTUM/ SI N( ALPHA) 

15.000 

12.000 


ALPHA 9.000 


6.000 

- 0.99999 - 

0.99999 



3.000- 

1.00000 - 1.00000 




0.000 

0.00000 





0.50000 0.75000 

1.00000 

1.50000 

2.00000 




H/0 


/-COMPONENT M0MENTUM/SIN1 ALPHA) 




15.000 



- 

0.99999 

12.000 



- 0.99999 - 

1.00000 

ALPHA 9.000 

- 

0.99999 

- 0.99999 - 

1.00000 

6.000 

- 0.96906 - 

1.16371 

- 0.99999 


3.300- 

0.7S15S - 1.06855 - 

0.99999 



0.000 

0.00000 





0.50000 0.75000 

1.00000 

1.50000 

2.00000 




H/D 


Z -COMPONENT MOMENTUM/SINt ALPHA) 




15.000 



- 

1.00657 

12.000 



- 1.01113 - 

1.00790 

ALPHA 9.000 

- 

0.87519 

- 1.01632 - 

1.01000 

6.000 

- 0.75791 - 

1.13051 

- 1.02156 


3.000- 

0.59022 - 0.99070 - 

1.07225 



0.000 

0.00000 

0.00000 




0.50000 0.75000 

1.00000 

1.50000 

2.00000 




H/D 


/-COMPONENT MOMENTUM/SINl ALPHA) 




15.000 



- 

1.00085 

12.000 



- 0.96033 - 

1.01273 

ALPHA 9.000 

- 

0.62151 

- 1.00971 - 

1.01775 

6.000 

- 0.55061 - 

0.99602 

- 1.03662 


3.000- 

0.39153 - 0.86778 - 

1.06332 



0.000 

0.00000 

0.00000 




0.50000 0.75000 

1.00000 

1.50000 

2.00000 




H/0 


Z-COMPONENT M0MENTUM/SIN1 ALPHA) 




15.000 



- 

0.87263 

12.000 



- 0.72169 - 

0.99729 

ALPHA 9.000 

- 

0.37006 

- 0.92051 - 

1.02867 

6.000 

- 0.36967 - 

0.78163 

- 1.03995 


3.000- 

0.18936 - 0.67969 - 

0.97979 



0.000 

0.00000 

0.00000 




0.50000 0.75000 

1.00000 

1.50000 

2.00000 




H/0 


Z -COMPONENT M0MENTUM/S1NI ALPHA) 




15.000 



- 

0.59667 

12.000 



- 0.66072 - 

0.88861 

ALPHA 9.000 

- 

0.15981 

- 0.73887 - 

1.01775 

6.000 

- 0.17319 - 

0.55339 

- 0.99001 


3.000 

0.00816 - 0.50772 - 

0.83373 



0.000 

0.00000 

0.00000 




0.50000 0.75000 

i. 00000 

1.50000 

2.00000 


H/0 

’• Normalized J values to show the effective development length 
z 

for the fully developed nozzle exit condition (J / J ( 0 ) sin a). 

z 
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Table 1. Radius r as a function of x/d to define the 0. 1 isotach 

(u (x, r) / u (o) = 0. 1) for the fully developed and uniform 
nozzle exit conditions. (The data shown are from those cases 
for which the measured x-component momentum flux was 
constant to within - 3%. ) 

Fully Developed 


a 

h/d 

x/d=l 

2 

3 

4 

5 

3 

0.75 

1.25 

1.41 

1.61 

1.82 

2.07 

3 

1 

1.21 

__ 

1.62 

__ 

2.01 

3 

1. 5 

1.25 

1.43 


1.85 

2.0 

9 

1 

1.26 

.. 

1.61 

.. 

__ 

12 

l 

1.21 

1.40 

1.63 



15 

1 

1.20 

1.38 

1.62 

.. 


15 

1. 5 

1.22 

1.42 

1.63 

1.83 


15 

2 

1.28 

.. 

1. 72 



Average 


1.23 

1.41 

1.63 

1.83 

2. 02 

Uniform 







0 

0. 75 

0. 9 

1. 03 

1.13 

1. 25 

1.40 

0 

1 

0.94 

1. 05 

1.20 

1.32 

.. 

3 

0. 75 

0.9 

1.01 

1.11 

1. 26 

1.40 

3 

1 

0.9 

1.0 

1.12 

1. 25 

1.38 

6 

1.5 

0.91 

1.03 

1.17 

1.30 

.. 

9 

1 

0.92 

1.03 

1. 16 

1.30 

_ 

15 

1 

0.93 

1.02 

.. 

.. 

«• _ 

Average 

« — 

0.91 

1.02 

1. 15 

1.28 

1.39 
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Table 2. Summary of percent widths, distance to maximum pressure 
and isobar which aligns with free isotach intersection. 


Width of 0. 1 x/d distance 

free isotach from x=o to 

intersection p max 

contour as a 

percentage of 

the zero isobar 

width recorded 

at the x- location 

of the maximum 

pressure 


p 

^agree^ 

X. p u(0) sin a 
all x I0~ c 


x/d distance 
from x=0 to 
0. 1 free 
isotach 
intersection 


a 

h/d 

u 

fd 

u 

fd 

u 

fd 

u 

fd 

3 

0. 5625 

.. 

74 

.. 

0.7 


0. 57 

*■» — 

0. 321 


0.667 

.. 

80 


2. 3 

.. 

0. 31 

w — 

1.003 


0.75 

0.7 5 

80 

4. 5 

2.9 

0. 35 

0. 266 

1.663 

1. 544 


1 

0.68 

78 

5.9 

5.7 

0.38 

0.266 

3. 543 

3. 176 

6 

0.75 

0.69 

76 

1.9 

2. 1 

0.36 

0. 30 

1. 16 

1. 18 


1 

0.75 

71 

4.7 

4. 2 

0.29 

0.30 

2. 5 

2. 324 


1. 5 

0.64 

68 

6.4 

7 

0.20 

0. 22 

5. 178 

4. 737 

9 

1 

0. 64 

mm 

3.2 

3.3 

0. 38 

0. 27 

1. 91 

1.81 


1. 5 

0. 65 

70 

5.7 

6.2 

0. 25 

0. 22 

3. 715 

3. 979 


2 


66 


7. 5 

.. 

0. 18 

. - 

5.619 

12 

1 

__ 


__ 

_ ^ 




- • 

• - 


1.5 

0. 58 

65 

5 

5. 1 

0. 28 

0. 224 

3. 202 

3. 029 


2 


65 


7.2 

.. 

0. 22 

• • 

4. 59 

15 

1 

0.46 

.. 

2 

— M 

1.0 


1. 247 



1.5 

__ 

.. 

.. 

m 


.. 

- - 

-- 


2 


61 


6 


0. 24 

4.06 

3. 855 

■ 

2 

0.3 


3.4 


0.65 

__ 

1. 978 

_ — 


3 

0.36 


5. 2 


0.61 


3.429 

- — 

60 

2 

0. 25 


1.9 


0.92 

_ _ 

0. 37 
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vO 


Table 4. Data from the Stagnation Point Investigation. Uniform Exit Condition v (4, 0, 0) = 


[p(x,0,0) - p atm ] 




/2 


P'd* 


X z — 

v ... (j"» 0» j - ) recorded from the hot-wire anenometer. 
w a a 


Aj of traverse and Aj to align v and v^ are indicated A^ = Aj (tan a + p)" 1 . 


a = 3 degrees, h/d = 1 


a = 1 

9 degrees, 

h/d = 1 


a =15 degrees 

, h/d = 

1 

Af =°. 

0328. A§ = 0 

.249 


4 - 

0.0395, A^ 

= 0. 164 


*5 = 

0.0 308, Ag 

• = 0.0871 

X 

d 

v p (fps) 

X 

d 

v w < f P 8 > ( 

X 

d 

Vp(fps) 

X 

d 

v w ( f P«) 

X 

d 

Vp(fps) 

X 

d 

vjfps) 

2.291 

0.000 

2. 166 

8.106 

0.875 

7.393 

0.958 

2.421 

0.000 

3.031 

0.000 

0.728 

2.458 

0.000 

2.333 

10.267 

1.04 1 

3.201 

1.125 

1.956 

0.166 

2.279 

0.166 

1.760 

2.625 

0.000 

2.500 

12.149 

1.208 

3.201 

1.291 

2.528 i 

0.333 

3.160 

0.333 

1.905 

2.791 

0.000 

2.666 

14.374 

1.374 

0.001 

1.458 

5.184 

0.500 

3.284 

0.500 

1.841 

2.958 

0.000 

2.833 

16.430 

1.541 

1.848 

1.624 

7.270 

0.666 

2.279 

0.666 

1.664 

3.124 

0.000 

2.999 

18.488 

1 .708 

4.802 

1.791 

10.428 

0.833 

3.845 

0.833 

0.606 

3.291 

0.000 

3.166 

20.618 

1.874 

10.496 

1.958 

14.903 

0.999 

3.322 

0.999 

4.838 

3.458 

0.000 

3.333 

22.374 

2.041 

13.862 

2.124 

19.776 

1.166 

13.290 

1.166 

10.633 

3.624 

0.000 

3.499 

24.077 

2.208 

17.534 

2.291 

24.703 

1 .333 

24.829 

1.333 

19.937 

3.791 

0.000 

3.666 

26.123 
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20.056 
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30.275 

1.499 

34.430 
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0.000 

3.833 
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2.541 
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2.624 

34.888 

1 .666 

40.723 

1.666 

40.350 

4.124 

0.000 

3.999 

29.605 

2.708 

24.711 

2.791 

39.815 

1.833 

44.38 5 

1.833 

50.807 

4.291 

0.000 

4.166 

31.387 

2.874 

25.577 

2.958 

43.902 

1 .999 

46.185 

1.999 

63.566 

4.458 

0.000 

4.333 

32.986 

3.04 1 

26.528 

3.124 

47.981 

2.166 

45.769 

2. 166 

81.836 

4.624 

0.000 

4.499 

35.254 

3.208 

27.352 

3.291 

51.562 

2.333 

45.031 

2.333 

93.062 

4.791 

0.000 

4.666 

36.443 

3.374 

27.305 

3.458 

54.864 

2.499 

43.768 

2.499 

97.307 

4.958 

0.000 

4.833 

38.389 

3.541 

27.321 

3.624 

58.891 

2.666 

42.056 

2.666 

100.325 

5.124 

0.000 

4.999 

39.859 

3.708 

26.816 

3.791 

60.690 

2.833 

40.283 

2.833 

103.503 

5.291 

0.000 

5.166 

41.628 

3.874 

26.302 

3.958 

62.490 

2.999 

38.768 

2.999 

104.968 

5.458 

0.000 

5.333 

42.329 

4.041 

26.188 

4. 124 

64.608 

3.166 

36.255 

3.166 

108.305 

5.624 

0.000 

5.499 

44.399 

4.208 

25.893 

4.291 

66.757 

3.333 

35.725 

3.333 

109.831 

5.791 

0.000 

5.666 

46.283 

4.374 

25.037 

4.458 

67.454 

3.499 

34.663 

3.499 

110.896 

5.95 8 

0.000 

5.833 

48.111 

4.541 

25.003 

4.624 

68.610 

3.666 

33.275 

3.666 

111.201 

6. 124 

0.000 

5.999 

50.033 

4.708 

24. 187 

4.791 

70.441 

3.833 

32.631 

3.833 

111.751 

6.458 

0.000 

6.333 

50.960 

5.041 

23.048 

5.124 

71.489 

4.166 

30.380 

4.166 

111.474 

6.791 

0.000 

6 • 666 

52.007 

5.374 

21.927 

5.458 

73.076 

4.499 

29.559 

4.499 

111.481 

7.124 

0.000 

6.999 

52.883 

5.708 

20.890 

5.791 

74.197 

4.833 

28.685 

4.833 

111.242 

7.458 

0.000 

7.333 

54.158 

6.041 

21.475 

6.124 

74.837 

5.166 

27.382 

5.166 

111.224 

7.791 

0.000 

7.666 

54.311 

6.374 

19.992 

6.458 

75.387 

5.499 

27.320 

5.499 

110.310 

8.124 

0.000 

7.999 

55.447 

6.708 

19.341 

6.791 

75.374 

5.833 

25.420 

5.833 

109.738 

8.45 8 

0.000 

8.333 

55 .b 1 1 

7.041 

18.713 

7. 124 

75.130 

6.166 

23. 184 

6.166 

108.238 
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Figure la. Externally blown flap STOL. 
aircraft configuration. 



Figure lb. The round- jet/plane- wall flow field, 
coordinate system and nomenclature. 



Figure 2. Detail of uniform flow nozzle 
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Figure 12. 
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a) x/d s» 0 


Figure 28. Isotach contours a = 3, h/d = 0. 75, uniform. 
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Figure 29. Isotach contours a = 9, h/d = 1, uniform. 
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d) x/d = 3 
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f) x/d = 5 
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Figure 30. 


Isotach contours a = 15, h/d = 1, uniform. 
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SIDE FRONT 

Figure 31. Control volume foi the analysis of the round- jet/ plane -wall 
flow field. 
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Figure 32. Normalized centerline pressure p/p\ u(o) sin a and 
the normalized curvature of the jet's momentum flux 
centerline Kd/sina . The x/d origin for each h/d value 
represents the intersection point of the 0. 1 free isotach. 
a r 3 degrees uniform exit condition 





0 1 2 3 x/d 4 5 6 7 


Figure 33. Normalized centerline pressure p/p X u(o)^ sin a and 
the normalized curvature of the jet's momentum flux 
centerline Kd/sina . The x/d = 0 origin for each h/d value 
represents the intersection point of the 0. I free isotach. 
a = 3 degrees fully developed exit condition 
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Normalized centerline pressure p/p X. u(o) sin a and 
the normalized curvature of the jet's momentum flux 
centerline Kd/sin a . The x/d = 0 origin for each h/d value 
represents the intersection point of the 0. 1 free isotach. 

a = 6 degrees fully developed exit condition 




p/p X u(o) si 
Kd/sin q 


Figaro 



36. 


Normalized centerline pressure p/p X u(o) sin a and 
the normalized curvature of the jet's momentum flux 
centerline Kd/sin a . The x/d = 0 origin for each h/d value 
represents the intersection point of the 0. 1 free isotach. 
a = 9 degrees uniform exit condition 


Kd/si 


si C 




h/d = 1.0 


Normalized centerline pressure p/p X u(o) 1 ' sin a and 
the normalized curvature of the jet's momentum flux 
centerline Kd/sin a . The x/d = 0 origin for each h/d value 
represents the intersection point of the 0. 1 free isotach. 

a - 9 degrees fully developed exit condition 



Kd/sin 



Figure 38. Normalized centerline pressure p/p \ u(o) sin a and 
the normalized curvature of the jet's momentum flux 
centerline Kd/sin a . The x/d = 0 origin for each h/d value 
represents the intersection point of the o. 1 free isotach. 
a =12 degrees uniform exit condition 
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p/p \ u(o) 
Kd/sin 


Figure 39. Normalized centerline pressure p/p \ u(o) 4 " sin a and 
the normalized curvature of the jet's momentum flux 
centerline Kd/sin a . The x/d = 0 origin for each h/d val u 
represents the intersection point of the 0. 1 free isotach. 

a =12 degrees fully developed exit condition 


Figure 40. Normalized centerline pressure p/p\ u(o) sin a and 
the normalized curvature of the jet's momentum flux 
centerline Kd/ sin a . The x/d = 0 origin for each h/d value 
represents the intersection point of the 0. 1 free isotach. 
a =15 degrees uniform exit condition 
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d = 2. 


Figure 41. Normalized centerline pressure p/p\ ufo)*" sin a and 
the normalized curvature of the jet's momentum flux 
centerline Kd/sin a . The x/d = 0 origin for each h/d 
value represents the intersection point of the 0. 1 free isotac 
a =15 degrees fully developed exit condition 


1 4Q 


Figure 42. Normalized centerline pressure p/p \ ufo) 4 " sin a and 
the normalized curvature of the jet's momentum flux 

centerline Kd/sin a . The x/d = 0 origin for each h/d value represent 
the intersection point of the 0. 1 free isotach. 

a = 30 degrees uniform exit condition 
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0 1 2 3 x/d 4 5 6 7 8 

Figure 44. Kd/sin a for the conditions a = 3 degrees, uniform, and the indicated h/d values 
The abscissa of the separate curves has been shifted such that the intersection 
of the 0. 1 free isotach occurs at the origin. 



values. The abscissa of the separate curves has been shifted such that the 


intersection of the 0. 1 free isotach occurs at the origin. 




o 1 2 3 x/d 4 5 6 7 8 

Figure 46. Kd/sin a for the conditions a = 3 degrees, fully developed, and the indicated h/d 

values. The abscissa of the separate curves has been shifted such that the 
intersection of the 0. 1 free isotach occurs at the origin. 
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values. The abscissa of the separate curves has been shifted such that the 
intersection of the 0. 1 free isotach occurs at the origin. 
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of the 0. 1 free isotach occurs at the origin. 
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intersection of the 0. 1 free isotach occurs at the origin. 
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Fig are 5 3. Kd/sin a for the conditions a =15 degrees, fully developed, and the indicated h/d 
values. The abscissa of the separate curves has been shifted such that the 
intersection of the 0. 1 f? isotach occurs at the origin. 
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Figure 54. Kd/sin a for the conditions a = 30 degrees, uniform, and the indicated h/d values. 

The abscissa of the separate curves has been shifted such that the intersection 
of the 0. 1 free isotach occurs at the origin. 






Figure 55. Kd/sin a for the conditions a = 60 degrees, uniform, and the indicated h/d values 


The abscissa of the separate curves has been shifted such that the intersection 
of the 0. 1 free isotach occurs at the origin. 


Axisymmetric region, azimuthal vorticity of the approaching jet flow. 

Buffer region, Streainwise vorticity production from reorientation of the 
azimuthal vorticity field 8u/39 ). 

Near wall region, streamwise vorticity from flux at wall ( y /9z) 
and vortex stretching (u 9 u/9 y) . 



Figure 56. Schematic representation of general isotach distribution 
showing the axisymmetric region in the undisturbed flow, 
the buffer region and the near wall region. 
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Normalized flux of 
streamwise vorticity 



Figure 57a. Flux of streamwise vorticity into the flow as a result of the 

surface pressure distribution, fully developed cases (a degrees, h/d). 
Note, 15, 1, u is for the uniform exit condition. 

NOTE: Abscissa for data shifted such that 0.1 free isotach aligns with the 
position x/d = 0. 


Normalized flux of 
streamwise vorticity 



Figure 57 b. Flux of streamwise vorticity into the flow as a result of the 

surface pressure distribution, fully developed cases (a degrees, h/ri). 
NOTE: Abscissa for data shifted such that 0. 1 free isotach aligns with 
the position x/d = 0. 
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Figure 59. Physical characteristics of oblique jet impingement. 
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Note: For the comparison of the hot-wire and pressure data, was 

shifted in x by an amount v w ( x . z) = v^fx+Ax, 0) where Ax = z w /tan a + |3. 


Figure 60. Experimental technique for the acquisition of the velocity and 
surface static pressure data. 
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Figure 61. Measured, v , and equivalent stagnation streamline velocities (as 

inferred from the surface pressure), v , to determine the stagnation 
point for the conditions a - 3 degrees, h/d = 1, uniform. 
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v = v measured at z/d = 0.0395 and 
A shifted by Ax = Az/tan (a + P) 
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Figure 62. Measured, v , and equivalent stagnation streamline velocities (as inferred 
from the surface pressure), v , to determine the stagnation point for the 
conditions a = 9 degrees, h/i- 1, uniform. 
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Figure 6 3 



Measured, v , and equivalent stagnation streamline velocities (as inferred 
from the surface pressure), v , to determine the stagnation point for the 
conditions a = 15 degrees, h/d 1, uniform. 
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Figure 64. Comparison cf for two z values, a 


15 degrees, h/d = 1, uniform. 



